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ABSTRACT 

We present ^^CO, ^^CO and C^^O (J=3-2) observations of a new cluster of 
outflows in the Vulpecula Rift with HARP-B on the JCMT. The mass associated 
with the outflows, measured using the ^^CO HARP-B observations and assuming a 
distance to the region of 2.3 kpc, is 129 M©, while the mass associated with the dense 
gas from C^^O observations is 458 Mq and the associated sub-millimeter core has 
a mass of 327 ±112 independently determined from Bolocam 1.1mm data. The 
outflow- to- core mass ratio is therefore ~0.4, making this region one of the most efficient 
observed thus far with more than an order of magnitude more mass in the outflow 
than would be expected based on previous results. The kinetic energy associated with 
the flows, 94x10^^ ergs, is enough to drive the turbulence in the local clump, and 
potentially unbind the local region altogether. The detection of SiO (J=8— 7) emission 
toward the outflows indicates that the flow is still active, and not simply a fossil flow. 
We also model the SEDs of the four YSOs associated with the molecular material, 
finding them all to be of mid to early B spectral type. The energetic nature of the 
outflows and significant reservoir of cold dust detected in the sub-mm suggest that 
these intermediate mass YSOs will continue to accrete and become massive, rather 
than reach the main sequence at their current mass. 

Key words: ISM: jets and outflows, stars: winds, outflows, molecular data, 
stars:formation, stars: pre-main-sequence 



1 INTRODUCTION 

Molecular outflows are observed towards stars of all masses 
I during their formation (jArce et al.ll2007l l, and a re generally 
, assoc i ated with both acti ve accretion (see e.g. Chu rchwell 

Il999l : IPudritz et aP l2007l ) and loss of angular momentum 
from the star /d i sk system tha t they originate from (e.g. 

iTomisakal |2000'; 'Ya mada etHI I2OO9I ). Early observations 
suggested that different scenarios for outflow generation 
might be required for the low and high mass regimes (see e.g. 

[Richer et al.ll2000 . and references therein.). However, higher- 
resolution studies of massive outflows have shown that these 
differences were exaggerated by the low-resolution single- 
dish observations used and that their properties transition 
relativ ely smoothly as a function of mass fe.g. lBeuther et aP 

l2002d li). 

In addition to being signposts of active star forma- 
tion, the mass, energy and angular momentum transported 
by molecular outflows can also have a significant impact 
on the surrounding molecular clouds in which young stars 
form. This can b e through the d r iving of turbulence within 
the cloud (e.g. iMatzned l2007l : ISrunt. Hever Mac Lowl 
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I2OO9I I. at least on small scales, or, for particularly power- 
ful outflows, even contribute to the dispersal of the cloud 
and removal of the res ervoir for further star formation 
(|Matzner McKedlioOol ) . A re cent study of the ou tflows in 
the Perseus molecular clou d by Curtis et al. (|201Q| . see also 
iHatchell k DunhamI 12009' and "A rce et al.ll2010h has shown 
that outflows can contribute a significant amount of en- 
ergy to driving turbulence and cloud dispersal for whole 
molecular clouds, though they may not always be ener- 
getic enough to dominate these processes on their own. 
Recent investigations into how the energy and momen- 
tum transported by outflows is translated to the cloud 
fe.g.lBaneriee Klessen FendtllioOTi : iNakamura Lill2007l : 
ICarroll et al.ll2009l ) show efficiencies which can be strongly 
scale- dependent, so even if outflows in a given cloud are not 
powerful enough to disrupt the whole cloud complex, they 
may be significant on the intermediate cluster scales. 

In this paper we present recent observations of a new 
outflow discovered at ^ = 59.6375°, b = -00.1875° in 
the ^^CO (J=l-0) observations of the Exeter FCRAO 
Five- College Radio Astr onomical Observatory) CO Survey 
Mottram k Bruntll20ld Brunt et al., 2011, in prep.), which 
have a spatial resolution of 45^^ and a velocity resolution of 
^0.15 kms~^. The systematic velocity of the cloud emission 
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Figure 1. Infrared three-colour GLIMPSE ("Beniamin et al.l l2003h IRAC 3.6 /im (blue), 4.5 idui (^reen) and 8.0 fim (red) image of part 
of the Vulpecula Rift with overlaid contours showing integrated Exeter-FCRAO Survey (|Mottram B runt 2010, Brunt et al., 2011, in 
prep.,) -"^-^CO emission. The contour levels for the integrated emission are 25, 50, 75 and 100 Kkms~^, and are coloured blue, green, 
orange and white respectively. The cyan dashed square indicates the region around the newly discovered outflow which will be explored 
in this paper, while the yellow dashed circle indicates the original identiflcation of the evolved H ll region Sh2-86 by Sharpless ( 1959). The 
magenta crosses indicate the positions of the OB stars identifled towards the cluster NGC6832 by Massev. Johns on &: Degioia-Eastwoodl 

1995), while the magenta dashed circle shows the approximate region containing the main cluster as identifled bv lBica, Bonatto Dutral 

2008l v 



associated with the outflow is ^27.5 kms ^, and the out- 
flow is esse ntially unresolved in these data. A three-colour 
GLIMPSE (jBeniamin et al.ll2003l l map, with contours show- 
ing the integrated FCRAO ^^CO (J=l-0) emission in the 
region and the outflow position marked is presented in Fig- 
ure [T] 

A sub-mm/ mm core ident if ied byTchapin et al.1 (|2008l l 
using BLAST (Pascale et all 12008) at ^ = 059.6331, 
b = —00.1906 (their V52), who derive a core mass of 
390 =b 75 M© and luminosity of 3850 =b 610 L©, is as- 
sociated with the outflow. The core is als o detected in 
the B olocam Galactic Plane Survey (BGPS ^Aguir re et al.l 
l201lh at 1.1mm which, using the median gas temperature 
of 27 K (obtained during outflow property calculation in 
section El , the source in tegrated flux and equation 10 from 
iRosolowskv et al.l(|2010h . gives a core mass of 327 ±112 M©. 

The outflow also lies near several mid/far-IR 
bright sources as shown in t he three-colour com bina- 
tion GLIMPSE, MIPSGAL (jCarev et al.l l2009l ) and 



Hi-GAL (i Molinari et~al1 l2010bl lah images in Figure El 
which are part of the Vulpecula Rift molecular cloud 
com plex. The region is near the large Hll region Sh- 
86 (IsharplessI Il959ll and the young clust er NGC 6832 
('Mas sev. Johnson Degioia-Eas twoodl Il995l . see Figure [T]), 
part of the VulOB 1 OB association. There have been some 
suggestions that the whole region is undergoing sequen- 
tial star f ormation (e .g . (Ehlerova. Palous Huchtmeierl 
l200lh . but iBillot et al.l (|2010l ) find no evidence for this 
in the YSO population as obtained from GLIMPSE and 
MIPSGAL observations. In t erms of the early phases of 
star formation and evolution, iBillot et al ] (|201lh found an 
increase in clustering of sources as wavelength decreases in 
the five Herschel bands (70 /im, 160 /xm, 250 /xm, 350 /xm, 
500 fim) of the Hi-GAL Science Demonst ration Phas e 
(SDP) observations of this region (Molinari et al.l l2010al ). 
While detailed identification of evolutionary phases has 
not been completed, this tentatively indicates that warmer 
sources are more likely to be in clusters within the region. 
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Figure 2. Infrared three-colour images of the region around the outflows. Left: GLIMPSE IRAC 3.6 /im (blue), 4.5 jum (green) a nd 
8.0 /im (red). Centre: GLIMPSE IRAC 8.0 /im (blue), MIPSGAL (Carey et al. 2009) 24 /im (green) and Hi-GAL (,Molinari et al.ll2010bl lah 
PACS 70 /im (red). Right: Hi-GAL PACS 70 /im (blue), PACS 160 /im (green) and SPIRE 250 /im (red). 



The molecular cloud properties of the whole Vulpecula 
Rift as observed in the Exeter FCRAO CO Survey will be 
discussed in an upcoming paper (Mottram et al, in prep.). 

The kinematic dis tance relating; to the systematic out- 
flow velocity, using the lBrand Blit d (|l993l ) rotation curve, 
is 2.7 kpc if the source is at the near distance and 5.9 kpc 
if at the far distance. The outflow is most likely as- 
sociated with the Vulpecula Rift molecular cloud com- 
plex, which infrare d extinction maps suggest is at the 
near distance (e.g. iRusseil et al.l l201lh . The open cluster 
NGC 6823 also lies within the cloud complex, for which 
iMassev. Johnson Degioia-EastwoodI (p. 995) derive a pho- 
tometric dista nce of 2.3 kpc. We ther efore follow the ap - 
proach of both lChapin et aP (|2008l ) and lRusseil et al.l 
in assuming a distance of 2.3 kpc to the region. 

It is perhaps surprising that an outflow associated with 
such a massive core has not been discovered sooner, but this 
might be because it only has an upp er limit at 25/im in the 
IRAS point source catalogue fPSC, [Beichman et al.|[l988l l 
which would have caused it to be excluded from the sele ction 
criteria of previous studies (e.g. Beuth er et aI1l2002d , who 
used the sample of Sridhar an et al. 20o3)"! 

We begin by describing the JCMT observations in sec- 
tion [21 after which we describe the method used to identify 
the outflow velocity windows and calculate outflow proper- 
ties in section [3l In section [H we present results and analysis 
of the outflow properties, as well as those relating to the lo- 
cal molecular material and the young stellar objects (YSOs) 
most probably associated with the outflow. Finally, we pro- 
vide a short discussion of the wider context of our results 
and reach our conclusions in sect. (5] 



2 JCMT OBSERVATIONS 

Mapping observations of ^^CO, ^^CO and C^^O (J=3-2) 
and single pointing grid observations of HCO^, H^^CO^ 
(J=4-3) and SiO (J=8-7) were obtained with HARP-B and 
the ACSIS autocorrelator (jBuckle et al.ll2009l l at the James 
Clerk Maxwell Telescope (JCMT), Mauna Kea, Hawaii on 
the V and 3^^ of November 2009 as service proposal 
S09BU01. The spatial resolution of the JCMT at the ob- 



served frequen cies is ^14^^ wi th a main beam efficiency 
r]MB = 0.61 (B uckle et aTll2009l l. Though HARP-B consists 
of a 4x4 pixel array of receivers, four of these pixels (HOG, 
HOI, H02 and H14) were not operational at the time of ob- 
serving. For each of the CO lines, four sets of observations 
in position-switching raster mode were undertaken with a 
quarter array shift (29.1^^) between each scan. A 90° scan 
direction change was also implemented between each pair of 
observations, after which the array was offset by an eighth 
(14.6^^) in both the x and y plane of the array and the process 
repeated. The resulting basket- woven maps cover an area of 
4^x4' centred on G059.6331-00.1906 (a = 19^43^49.68", 
6 = +23°28'38.3'' and Ylsr = +25 kms"^). The correlator 
settings and observing times are summarised in Table [1] 

The data were initially reduced using the Starlink pro- 
gram ORAC-DR, however the broad wings of the outflow were 
within the baseline region determined by this routine. The 
data were therefore rebaselined and smoothed with a gaus- 
sian kernel to lower velocity resolution in order to reduce 
noise using python scripts. During the ^^CO and C^^O ob- 
servations, one of the HARP-B detectors (H04) had a partic- 
ularly bad baseline, so after examination these data were re- 
reduced with this detector excluded. All observations were 
then converted to the main beam temperature scale (Tmb). 



DETERMINATION OF OUTFLOW 
PROPERTIES 



Previous studies of outflows (e.g, 



iBeuther et aD l2002d : IStoiimirovic et al.l 



Bally et al.l Il999l : 

y) have tended 



200i 



to define global velocity windows, often by eye, to apply to 
all spectra containing outflow emission, or to a single spa- 
tially integrated spectrum. While this may work well for 
simple regions with only one outflow, using a single velocity 
window for all spectra across a region assumes that the cloud 
systematic velocity and line- width do not vary with position 
and that the outflow has a roughly similar velocity extent in 
all spectra. In the case of our observations this is not a safe 
assumption, as the systematic velocity and the full- width 
half-maximum (FWHM) both vary by '-^2kms~^ across the 
outflow region and the velocity extent varies strongly with 
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Table 1. Observational Setup 



Transition 


Rest frequency 
(GHz) 


Bandwidth 
(MHz) 


No. of channels 


Vel. Resolution 
(kms-i) 


Obs. Type 


On source time 
(min) 


CO (J=i-o) 


345.80 


250 


4096 


0.053 


Raster 


52 




345.80 


1000 


1024 


0.847 


Raster 


52 


i^CO (J=l-0) 


330.59 


250 


4096 


0.055 


Raster 


71 


C^^O (j=l-0) 


329.33 


250 


4096 


0.056 


Raster 


71 


HCO+ (J=4-3) 


356.73 


250 


4096 


0.051 


Grid 


10 




356.73 


1000 


1024 


0.821 


Grid 


10 


SiO (J=8-7) 


347.33 


250 


4096 


0.053 


Grid 


20 


H13CO+ (J=4-3) 


347.00 


250 


4096 


0.053 


Grid 


20 



position. We therefore set out below the method and cri- 
teria we have used to determine the velocity window and 
mass individually for each spaxel within which emission is 
attributed to the outflow. 

(i) Some ^^CO spectra show additional cloud emission at 
different systematic velocities to the main feature associated 
with the outflow (e.g. see Figure [3|). Gaussian fits were used 
to remove these, where possible, from the spectra before 
measurements of the outflow properties were performed. 
Strict acceptance criteria were used for the fit results so 
that the peak velocity and line- width did not vary too much 
from the expected values in order to ensure that these fits 
did not remove outflow emission. 

(ii) The maximum temperatures of both the ^^CO and 
^^CO spectra were measured using gaussian fits to the 
central cloud emission. 

(iii) The excitation temperature (Tex) was calculated for 
the spectrum, assuming that the ^^CO emission is optically 
thick and that the gas along a given line of sight can be 
characterised by a single excitation temperature. 

(iv) The ^^00 optical depth (r^^) was calculated, 
assumme: that the ^^00 emission is optically thin. 

(v) The systematic velocity of the cloud for the spectrum 
was taken to be the velocity of the peak of the gaussian 
fit to the ^^CO emission line, as the ^^CO spectra do not 
include a significant contribution from the outflow. 

(vi) The emission at velocities within 3cr of the peak of 
the gaussian fit to the ^^CO spectrum was associated with 
the cloud, rather than the outflow, in the ^^CO spectrum. 

(vii) We estimate the ^^CO to ^^CO isotopic ratio as 
a function of velocity (Ri2/i3('^)) for the spectrum, using 
a quadratic fit to the velocity window identified as cloud 
emission, in order to account for optical depth effects 
in the ^^CO spectr um. The maximum val ue of the ratio 
was limited to 62 (|Langer &; Penziad llQQSi ). the assumed 
intrinsic ratio when both transitions are optically thin. 

(viii) The ^^CO column density per velocity channel 
{N{^^C>0){v)) was calculated from the ^^CO spectrum, 
using Ri2/i3('^^) and r^^ to account for optical depth effects. 




Aw (km s"^) 



Figure 3. Top: Example ^^CO (black) and -"^^CO (green) spectra. 
The yellow line indicates a gaussian fit to a non-outflow emis- 
sion feature, which was removed when calculating the outflow 
properties. The inner and maximum outflow velocities are shown 
by the dashed and dash-dotted lines respectively, with the -"^^CO 
spectrum between these velocities shaded the appropriate outflow 
wing colour. Bottom: Gumulative integral of the mass in the out- 
flow wings versus velocity offset for the spectrum shown above. 
The error bars indicate the random error associated with the cu- 
mulative integrated mass. The total integrated mass, obtained as 
discussed in section [3] step |(xi)| is indicated by the black line, 
while the velocity at which the cumulative integrated mass flrst 
becomes greater than or equal to the total mass, used as the 
maximum outflow velocity, is marked by the dot-dashed lines. 
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(ix) The mass per velocity channel {mH2{v)) was calcu- 
lated from the ^^CO column density per velocity channel, as- 
suming that the ratio of the H2 column density to the ^^CO 
column density is 1.2x10^ (^Frerkin g. Lan ger WilsonI 
Il982l ). that the ratio of the ^^C O to ^^CO column density 
is 62 ([Langer Penzias Il993l), that the mean molecular 
weight of the gas is 2.8 ( Kauffmann et allfioQS ). and using 
data with a pixel size of 7.3" . 

(x) The cumulative mass per velocity channel was 
calculated separately for the red and blue outflow lobes as 
a function of the velocity offset (Av) from the systematic 
velocity of the cloud, starting at the inner outflow velocity 
defined in step (vi) and stopping at a velocity offset of 
57.5km s~^, determined to be free of outflow emission in all 
spectra (see Figure [3] for an example). 

(xi) We calculate the total integrated mass in the outflow 
lobe for the spectrum using a constant least-squares fit 
to the velocities larger than the point at which the ^^CO 
spectrum first becomes negative, i.e. has roughly reached 
the background level. The maximum velocity of the outflow 
lobe for the spectrum is then defined as the velocity at 
which the integrated mass first becomes greater than or 
equal to the total integrated mass (see Figure [3] for an 
example) . 

(xii) The momentum per velocity channel is calculated 
by multiplying the mass per velocity channel calculated 



in step (ix) by the velocity offset of each channel from 
the systematic cloud velocity and the kinetic energy per 
velocity channel by multiplying by half the velocity offset 
squared. The total momentum and kinetic energy in each 
lobe for the spectrum are then calculated as the integrals 
of the momentum and kinetic energy per velocity channel 
between the inner and maximum outflow velocities, as 



defined in steps (vi) and (xi 



While removing other cloud emission features in step (i) 
may have resulted in removal of a small amount of outflow 
emission, we consider such a conservative approach to be 
preferable to not removing this cloud emission. Examples 
of the regions where such emission has been identified for 
removal are indicated using yellow contours in the position- 
velocity plots in Figure SI while the outflow velocity win- 
dows are shown using red and blue contours respectively. 
The equations us ed for steps (iii)|(iv)|(viii) and (ix) can be 
found in lWilsonI (|2009l l. 

While the random uncertainty in the mass due to ther- 
mal noise is relatively small (e.g. see the error bars in Fig- 
ure [3j_2_th^se_aTeundoubtabty dominated by systematic er- 
rors. i Buckle et al ] ([2009) estimate the systematic calibra- 
tion errors for the JCMT to be ^20% and the various abun- 
dance ratios used are not particularly well know. In addition, 
the values used here were primarily measured in average 
ISM/molecular cloud conditions, so may not be representa- 
tive of the abundance in star formation regions where deple- 
tion, photodissociation and fractionation processes can vary 
depending on position, species and even isotope. Masking 
within 3(7 in step |(vi)| is reasonably conservative, as shown 
in Figure (4] and probably leads to some outflow emission at 
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Figure 4. Outflows identified from HARP-B data. Top left: In- 
tegrated 13 CO emission shown in greyscale with red and blue 
contours showing outflow emission detected in ^^CO for the same 
region as in Figure [2] The red and blue contour levels are 25, 
50, 75, 100 and 125 Kkms"-*^. Right, top and bottom: Position- 
velocity (PV) cuts through the -"^^CO data, indicated by diagonal 
numbered lines in the top right plot, shown in a square-root scale 
to emphasise low-scale emission. The levels of the black contours 
to the I^CO data are 0.8, 2.15, 5.6, 10.5, 17 and 25kms-\ while 
the yellow contours follow the same scale and indicate the regions 
of non-outflow emission removed using gaussian fits as discussed 
in the text. The red and blue contours indicate the red and blue 
outflow velocity windows. Bottom Left: Example -"^^CO (black) 
and -"^^CO (green) spectra for the pixel at the centre of the two 
PV cuts, with all lines and shaded regions having the same mean- 
ing as the top plot in Figure [S] 



low velocity offsets being removed, perhaps resulting in as 
large as a factor of 3 un derestimation of the outflow mass 
(jPownes &: CabritI bOC? !) . However, this approach is more 
accurate than either not correcting for ambient material, 
which can result in a significant overestimate, or attempt- 
ing a more complex correction without sufficient information 
(|Cabrit k Bertoutlll99Q ). Despite our conservative masking 
of the cloud emission, it is also possible that some of the 
emission considered as associated with the outflow is actu- 
ally due to local turbulent material. While it is difficult to 
estimate a numerical error related to these factors, by us- 
ing abundance ratios which are commonly used by other 
authors in the literature, we expect that we are not signif- 
icantly more susceptible than other studies. The properties 
presented below are probably accurate to a factor of a few. 
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Figure 5. HARP-B ■'^^CO channel maps, integrated in 6 kms~^ slices except for the four slices at near the systematic velocity, which 
are 2.5 kms"-*^ in order to better reveal the exclusion of systematic emission from the outflow. The region to considered when calculating 
the total mass in the region is indicated by the green contour, while the velocity integration of each slice is given in the top right corner 
of each plot. Those slices where at least 1% of the pixels lie within the outflow velocity window are denoted with a star and the text 
is colour-coded to the outflow wing. Contours of the total integrated emission in each outflow lobe are overlaid on the central velocity 
slices using the same contour levels as in Figure [3] The beam-size is shown in the bottom left of the plot. 
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4 RESULTS 

4.1 Outflow properties 

We performed the steps outlined in section [3] on all HARP 
^^CO spectra within a spatial radius of 10 pixels (^O.Slpc) 
of the centre of the source, as this region encompasses the 
whole outflow. A series of channel maps with the channels 
identified as contributing to the outflow indicated are shown 
in Figure [3 with the region where our automated outflow 
identification was undertaken shown by the green contour. 

It is clear from Figures [H and [5] that the single outflow 
originally identified in the FCRAO ^^CO (J=l-0) data is 
at least partially resolved into multiple components at the 
resolution of the HARP-B observations. By examining both 
the channel and integrated maps, the blue wing emission 
breaks up into 4 or 5 components, while 3 or 4 separate 
components are visible in the red wing, though it is diffi- 
cult to unambiguously identify the individual flows. In the 
following analysis, we tentatively identify four outflow com- 
ponents from the integrated emission maps in both the red 
and blue wings, as indicated in Figure (6] for which we have 
measured the mass, maximum relative velocity, momentum 
and kinetic energy, presented in Tabled In regions of over- 
lap, we assign spaxels in a mutually exclusive way. In order 
that our determination of individual flows does not overly 
affect our results, we also calculate the total properties of 
all the outflow emission. It is difficult to accurately estimate 
the angle of inclination between the line of sight and the 
outflow axis simply from our molecular observations, as the 
individual components are not well resolved or elongated, so 
we do not undertake an inclination correction for our veloci- 
ties, momenta and energies. Thus the values of the quantities 
calculated below are lower limits to the true values and we 
are unable to calculate dynamical timescales. 



4.2 The Ambient Molecular Cloud 

The molecular clump associated with the outflow lies within 
a radius of ~0.81 pc (10 pixel), shown by the green contour 
in Figure El and has a mass, measured from the ^^CO obser- 
vations, of 3655 M0, with 458 in dense gas as measured 
using the HARP-B C^^O (J=3-2) data and a ratio of the 
C^^O a nd ^^CO column densities (i.e. Nj^ ^CO) / N{C^^O)) 
of 500 (jPrerking, Langer & WilsonlflQSi ). The same veloc- 
ity window was used as for the outflow mass measure- 
ments, but without masking the line centre. Assuming that 
clump is spherical with a radius of 0.81 pc, we calculate 
its gravitational binding energy from the ^"^CO mass to 
be 1.4x10^^ ergs. We also calculate the turbulent energy 
(Eturb) of the clump using: 



{3 /16ln2)Mcioud AVfwum 



(1) 



from lArce GoodmanI (|200lh to be 1.7x10 ergs, 
where Vfwhm is the average ^^CO FWHM of 2.4 kms"^ 
and Mcioud is the ^^CO clump mass. The combined energy of 
the outflow, even without inclination correction, is therefore 
approximately enough to drive the turbulence in the clump, 
and would provide enough energy to unbind the clump if 
i<15°. 
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Figure 6. Red and blue outflow components indicated with 
coloured ellipses overlaid on the integrated -"^^CO emission map 
corresponding to that line wing, as shown with contours in Fig- 
ure [5] The numbers in the top right corner of each plot are the 
component numbers, also used in Table [2] The green contour is 
the same as in Figure [5] The four sources (sl-s4) associated with 
the clump containing the outflows are indicated in magenta. See 
sect. 14.31 for more details. 



4.3 Associated YSOs 

In order to explore the (proto) -stellar sources which are 
potentially associated with the detected outflows, we first 
note that four sources can be seen in Figure [2] which are 
probably associated with the molecular clump. Photome- 
try data for these sources was gathered from t he UKIDSS 
Galactic Plane Survey at J, H and K bands (Lucas e t al.l 
2008), the GLIMPSE s urvey in the four Spitzer IRAC bands 
|Beniamin et al.ll2003l ) and prelimin ary PACS 70 /xni and 
160 Mm from t he Hi-GAL survey (|Molinari et al.l l2010al : 
lElia et all l2010h . For the two sources not detected at all 
four IRAC bands, and for all sources from MIPSGAL 24 fim 
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Table 2. Measured outflow properties 



Component 


M 




Av 


max 




P 


E 






(Mo) 


(km 


s-1) 


(Mo 


km s"-*^) 


(10^5 ergs) 




r 


b 


r 


b 


r 


b 


r 


b 


1 


8.5 


30.3 


38.6 


38.9 


48.9 


196.3 


3.29 


15.93 


2 


9.7 


26.0 


50.5 


41.6 


74.9 


142.1 


8.94 


9.50 


3 


18.7 


5.8 


48.5 


42.6 


173.3 


44.8 


28.13 


4.51 


4 


22.9 


6.9 


46.9 


37.4 


175.7 


52.9 


19.02 


4.90 


total 


59.8 


69.1 






472.9 


436.1 


59.38 


34.84 



Table 3. Properties derived from SED fltting 



Source 




logio(L/LQ) 


si 


3.2 ± 0.9 


2.7 ± 0.2 


s2 


33.8 ± UA 


3.8 ± 0.2 


s3 


50.8 ± 9.4 


3.9 ± 0.2 


s4 


4.0 ± 0.7 


2.8 ± 0.2 


total 


91.8 ± 17.2 


4.2 ± 0.4 



images (|Carev et al.ll20o"9l ) aperture fitting photome t ry wa s 
performed to obtain flux data as in iMottram et"all |2010). 
Sources 2 and 3 are mildly saturated in the 24 /xm images, 
so the fluxes for these sources are treated as lower limits. In 
the Hi-GAL SPIRE 250 /xm, 350 /xm and 500 /xm images, 
the sources are too close together for reliable photometry 
to be obtained for the individual components. We therefore 
performed one aperture fitting photometry measurement en- 
compassing all the sources, and use the flux obtained as an 
upper limit for sources 2,3 and 4. Source 1 is not included 
in this process because it is not detected at 160 /xm. 

The spectral energy distributions (SEDs) of 
these sources w e re th en fit using the model fitter of 
iRobitaille etaP (12007 



^ with similar input parame- 
ters to those used in i Mottram et al.1 (HonB). We use 
a distance of 2.3=t0.5 kpc for all sources, where the 
distance error is set conservatively based on the differ- 
ence between the photometric and kinematic distances 
rather than the error on the photom e tric d istance from 
iMassev, Johnson Degioia-EastwoodI (|l995l ), which is 
much smaller. The SEDs and fits are show in Figure [3 
while the results are presented in Table [3] 

Despite the increase in resolution between the FCRAO 
and HARP-B data, it is still not simple to associate the 
sources directly with individual outfiow components. How- 
ever, sources 1, 3 and 4 all seem spatially related to the ob- 
served outfiows, partic ularly sources 3 and 4. B ased on val- 
ues given m Table 1 of iMottram et al.l (|2011al l. the sources 
have luminosities consistent with spectral types in the mid to 
early B type, with the most massive (S3) being a Bl— B0.5. 
The absence of detected radio continuum emi s sion in Red 
MSX Source (RMS) Survey (lUrauhart etlH I2OO8I] VLA 
5 GHz observations of this region (^Urquhart et al. I I2OO9I ) 
suggests that these sources are not powering Hll regions, 
so have yet to reach the main sequence. 



4.4 SiO, HCO+ & H^^CO+ 

The shock front and heating caused by molecular out- 
fiows interacting with surrounding cloud material strongly 
effects the chemistry of the surrounding region by re- 
leasing depleted molec ules like SiO from grain mantles 
into the gas phase (e.g. ICaselli, Hartquist Havneslll997l : 
ISchilke et al.l Il997h , increasing it's relative abundance by 
fact ors of ^10^ — 10^ with respect to the ambient material 
fe.g. lMartin-Pintado, Bachiller Fuentelll992l : iGarav et all 
Il998h . Given that the depletion timescale for SiO onto 
dust grains is relatively short in nor mal molecular cloud 
conditions , of order 10^-10"^ yrs (Mikami et aD Il992l : 
[Martin-Pin tado, Bachiller Fuentel ll992), the presence of 
this molecule is therefore a good indicator that an out- 
fiow is currently active, rather than being a fossil fiow. 
Our SiO observations are shown in Figure (8] where we de- 
tect broad emission (FWHM?^20 kms"^) towards the cen- 
tre of the outfiow, indicating that the outfiows in this region 
are indeed active and not left over from previous activity. 
The amplitude ob served is consiste nt with other observa- 
tions (e.g. Klaasse ^^Wilsonl l2007l ). The H^^CO+ obser- 
vations have a line-centre of 27.8 kms~^ and a FWHM of 
3.3 kms~^, consistent with the C^^O lines for the same po- 
sition, so the HCO+ and H^^CO+ emission is most likely 
associated with the core/envelope rather than the outfiow. 
The HC0+/H^3C0+ ratio for this location is in the range 
7—28, much lower than the expected abundance ratio for 
i2q^i3q q£ g2 indicating that HCO+ is optically thick and 
thus the double-peaked emission profile is almost certainly 
due to self-absorption. 



5 DISCUSSION & CONCLUSIONS 

The total mass within the outfiow components, as obtained 
from the HARP-B ^^CO (J=3-2) observations is 129 Mo, 
while the core mass derived from the BGPS 1.1mm dust 
continuum source is 327 ±112 Mo- The mass in dense gas, 
as measured from C^^O (J=3— 2) data is 458 Mo, which 
agrees reasonably well with the core mass. 

Using the BGPS core mass and the r elationship be- 
tween outfiow and core ma ss obtained b y iBeuther et al.l 
('2002c", their Fig 7.) to the Beuther et all (|2002bl ) sample 
results in an expected outfiow mass of ~10.3 Mo, an or- 
der of magnitude smaller than detected. The outfiow-to-core 
mass ratio is . 4, high er than any of the sources reported by 
iBeuther et al.l (|2002bl H) , though the total luminosity of the 
associated YSOs is not particularly high compared to their 
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A (Aim) 

Figure 7. SEDs of the four YSOs associated with the outflow 
material, located as indicated in Figure [6] Detections are indi- 
cated with fllled black circled while upper and lower limits are 
indicated with upwards and downwards pointing fllled black tri- 
angles respectively. The best flt model is shown as a black line 
while all other flts within a given Ax^ of the best flt are shown 
by grey lines. 



sample. This suggests that these new outflows are either 
more eflficie nt at entraining materi al or more energetic than 
those in the'Beuther et al.1 (|2002b) sample. If we assume an 
age for the outflows of 10^ yrs (5x10^ yrs), the mechan- 
ical luminosity of the combined outflows is 79L0 (2L0), 
resulting in a ratio of the mechanical to YSO luminosity 
of 0.5% (0.01%). In order to compare with the results of 
iBeuther etHI (|2002ch we must first account for the increase 




V (km s-i) 



Figure 8. HCO+ (black), R^^CO^ (blue) and SiO (red) spectra 
smoothed to 0.5 kms"-*^, 0.5 kms"-*^ and 2.5 kms"-*^ respectively 
for the central position of the outflow. The H^^CO"^ spectra have 
been multiplied by 5 and the SiO by a factor of 25 for improved 
visibility 



which results from using the maximum outflow velocity and 
total mass rather than the channel velocity and mass when 
calculating the kinetic energy (approximately a factor of 25). 
Once this is taken into account, the highest mechanical to 
YSO luminosity ratio in their sample is 0.07% for IRAS 
19410+2336 (G059.7833+00.0647 visible in the top left of 
Figure [1]) and the median ratio is 0.002%. 

There are several factors which could result in such high 
out flow- to- mass and mechanical-to-source luminosity ratios 
compared with previous studies. Firstly, it is possible that 
the additional driving source(s) for the outflows are heav- 
ily embedded within dense dust and gas core, though if this 
were the case it is surprising that there is not even a sug- 
gestion of this in Figure [2] as these data cover 3— 250/im. It 
could also be the case that despite our best efforts to be con- 
servative, some cloud emission is still being included within 
the outflow velocity windows, but our mass measurements 
would have to drop by an order of magnitude to fall in line 
with other sources. The gas might have a higher than normal 
CO to H2 ratio, leading to an overestimation of gas, though 
one has to ask why this particular region would be special 
in this regard. 

Finally, it could be that outflow activity from 
young massive stars is variable due to variable ac- 
cretio n, in a similar way to the senario suggested 
by Baraffe, Chabrier &: Gallardol (|2009l , see also 
'Baraf fe k Chabried 12010^ for low- mass stars, where 
sources spend ~1% of their time in a high- accretion 
(lO~^M0yr~^) phase and the rest accreting at a much 
lower rate (IO^^Mq yr~^). If the outflow from such a source 
was observed during or soon after the high accretion phase, 
a much higher mechanical luminosity would be observed 
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than after a long low-activity phase but the luminosity and 
core mass of the source would be similar. 

In summary, we have confirmed an outflow with mul- 
tiple components, probably associated with a group of four 
YSOs, in the Vulpecula Rift which has an outflow mass of 
129 M0 and a core mass of 327 =b 112 M©. The combined 
kinetic energy in the outflows (94x10^^ ergs) is enough to 
drive turbulence in the local clump and potentially unbind 
the region, depending on inclination and the efficiency with 
which that energy is translated to the cloud. The outflows 
are not only energetic but also certainly active, due to the 
SiO detection, and thus presumably the related mid to early 
B type YSOs are currently still undergoing a phase of major 
accretion. Given the large reservoirs of cold dust detected in 
the sub-mm towards these sources, most likely associated 
with dense envelopes, this phase will probably continue for 
some time. The relatively red SEDs and the energetic nature 
of the outflows point to the relative youth of these sources, 
so it is likely these YSOs are destined to become more mas- 
sive, rather than reach the main-sequence at their current 
mass. 

Overall, some variation in out flow- to- core mass ratio is 
to be expected due to variations in cloud conditions, accre- 
tion rates, entrainment efficiency and the powering source. 
However, the outflow mass, momentum and energy must 
be in some way related to the mass and age of the central 
source, and thus the energy able to be injected, as found em- 
pirically by iBeuther et al.1 (|2002d ). To find a young region 
where the outflows appear to be an order of magnitude more 
efficient than expected warrants further detailed investiga- 
tion, in order to explore how outflows entrain mass from and 
impart energy and momentum to the surrounding molecular 
material. We therefore intend to undertake further interfer- 
ometric observations at higher spatial resolution to study 
this region. 
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